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We discuss low-temperature properties of magnets with integer spin and large single- ion easy-plane 
anisotropy D in transverse magnetic field h. Considering the exchange interaction between spins as 
a perturbation and using the diagram technique we derive at h ~ D in the first nonvanishing orders 
of the perturbation theory thermal corrections to the elementary excitation spectrum, magnetization 
and specific heat. An expression for the boundary h c i(T) is found in the h-T plane between the 
I paramagnetic phase and that with the long range magnetic order. The effective interaction is 

derived between bosons near the quantum critical point h c i(0). The proposed theory describes well 
q experimental data obtained in NiCl 2 -4SC(NH 2 )2 (DTN). 
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Introduction. The topic of quantum criticality has received much attention in recent two decades. Of particular 
interest are quantum critical points (QCPs) which can be reached in experiments by varying easily controllable 
parameters such as external magnetic field, pressure, level of doping, etc. The equivalence between a spin system 
and a diluted gas of bosonic particles proved to be very useful in describing field-induced QCPs in magnets.^ This 
equivalence is revealed and exploited using appropriate representation of spin operators via bosonic ones. 

We discuss in the present paper properties of a system on a 3D lattice with an integer spin and large single-ion 
anisotropy which is described by the Hamiltonian 

n = dY,(sd 2 + \ E J iJ s i s j + hY, s *> (!) 

where D > is assumed to be much larger than exchange constants (D 3> J) so that the ground state at h = 
is paramagnetic (all spins are approximately in the state with S z =0). This system has at least two field- induced 
QCPs corresponding to transitions from the paramagnetic (at h — h c ±(T — 0)) and from the fully polarized (at 
h = h C 2(T = 0)) phases to other phases which nature depends on the detail of the exchange coupling and the lattice 
geometry. We propose in our recent paper 2 a bosonic integer spin representation that is convenient for discussion 
of the system ([I]) properties in the paramagnetic phase. Using this representation and considering the exchange 
interaction as a perturbation we find in RefP the spectrum of the Hamiltonian ([!]) in the paramagnetic phase at 
h = in the third order in the perturbation theory (hereafter referred to as expansion in terms of J/D for shot). We 
continue our study of the model (JlJ in the present paper and address its low-temperature properties in the vicinity 
of the QCP h — h c i(T — 0) using the proposed bosonic spin representation. Expressions are derived below in the 
first nonvanishing orders in J/D for thermal corrections to the elementary excitation spectrum, magnetization and 
specific heat. An expression is found for h c i(T) that is the boundary of the paramagnetic phase in the h-T plane. The 
effective interaction is derived between bosons near the QCP which can be observed experimentally. We demonstrate 
that the proposed theory describes well corresponding experimental data obtained in NiCl2-4SC(NH2)2 (DTNjP^ 
which is the most extensively studied compound having the paramagnetic ground state at h = and which is modeled 
by the Hamiltonian (JlJ . 

The magnetic subsystem of DTN consists of Ni ions with 5=1 and the Lande factor g = 2.26. Magnetic ions form 
a body-centered tetragonal lattice which can be viewed as two interpenetrating tetragonal sublattices. The exchange 
interaction between spins inside one sublattice is antiferromagnetic and strongly anisotropic: the exchange constant 
along the tetragonal hard axis (z axis) is much larger than those along x and y axes. Then, DTN is a quasi-lD 
material having two QCPs alP^ 

gfi B h° TN (T = 0) = 2.05T, (2) 
gfi B hg TN {T = 0) = 12.175 T (3) 

with a canted antiferromagnetic phase between them. It was found™ that the QCP h = h cl (0) belongs to the 3D 
BEC universality class: h c \(T) — h c i(0) oc T a with a s» 1.5. The strength of the effective interaction between long- 
wavelength boson was extracted in RefP from the measurements of magnetization and h c \{T). The specific heat at 
small T was measured in RefsP^l at h ~ h c \{0) and h ~ h C 2(0). 

We have shown in our previous papeil^ that Hamiltonian ([l]) describes well the experimentally obtained spectrum 
of elementary excitations in DTN at h = with the following set of parameters which differs from the conventional 
one (see discussion in RefP): 

D = 7.72 K, 

J z = 1.86 K, (4) 
J xy = 0.2 K, 
V = 0.1 K, 

where J z is the exchange constant along the chains, J xy is the exchange coupling constants between chains inside one 
tetragonal sublattice and V is the exchange coupling constant between neighboring spins from different tetragonal 
sublattices (which is proposed in RefP^). We use these parameters below for the experimental data analysis. 

Method and technique. It is convenient to use the following spin representation of integer S which is proposed in 
our previous studyP 



S? = b\bt - oJo(, (5) 
S+ = Sf 




c\ - c 2 a\a t ) a { (6) 




(S-frfoxS + i + bjfr.) / (g- Q t gi )(5 + i + a t ai ) 
1 + 0% V 1 + a)ai 
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where aj and bi are bosonic operators, ci = y/S(S + 1) and c 2 = \/S(S + 1) — ^/(S 1 — 1)(5 + 2)/2 > 0. Representation 
([5])-([6]) reproduces the spin commutation relations on the physical subspace which is constrained by the following 
additional term in the Hamiltonian (see discussion in RefP): 



i 

Substituting Eqs. ([5])-([6| into Eq. ([I]) and taking into account Eq. ^ one obtains for the Hamiltonian 

2 

^ = E [ £ ia(p) a p a p + £ ib(p) fo p 6 p] + E c ~\ 3 v ( fl P 6 -P + a P fe -p) 
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Pl+P2+P3=P4 

where J p = J2i Jije ipRii , N is the number of unit cells, 



(a\ ala 3 a,4 + b\ + [U — J3-1] a{&|a 3 & 4 
Ji (b\a\a\a4 + a\b\b\bi + a\a 3 a 2 bi + blb^a^J , 



eio,fc(P) = -° + y^P - 
c 2 

eu,fc(p) = £> + y^p + ^ 



(7) 

(8a) 
(8b) 
(8c) 

(9) 
(10) 



are spectra of a and b particles in the first order in J/D (here and below the number in the lower index of the spectrum 
indicates its order in J/D). It is convenient for the following to introduce Green's functions 



G a ,h(p) = -*( a p«p)> 

G b ,h(p) = -i{ b P b l)i 
F h (p) = -i(bl p al) 7 



(11a) 
(lib) 
(11c) 



where p = (uj, p) and a„ is the Fourier transform of a p (r). Dyson equations for one couple of these Green's functions 
have the form 



G a ,h(p) = G 0a ( P ) [1 + Za,h(p)GaAp) + ^h{p)F h {p)] 

F h ( P ) = Go b {-p) pfc(p)G , h (p) + E b , h (-p)F h (pj\ . 



(12) 



where Go a (p) = (w — ei a ,ft(p) , Gob(p) = (<*; — £ib,h{p) , S and II are normal and anomalous self-energy 

parts, respectively. Solving Eqs. (12) and the couple of equations for Gb,h{p) and F h {p) one obtains 



Ga,h(p) 

Gb,h(p) 
Fh{p) 



uj + E^hj-p) 
V{p) ' 

g + £q,ft(-.P) 

n h (p) 

2>(p) ' 



2?(p) = (w - eia(p) - E„, h (p)) (w + eib(p) + E 6lh (-p)) + \H- h {p)\ 2 . 
Spectra of a and b particles are given by equations 

■D(e a>h (p), p) = 0, £>(-e 6 ,h(p), p) = 0. 
T = artrf /i = 0. Because a and b particles are equivalent at h = we have 

G a ,h=o{p) = G b ,h=o(p) = G{p), 

n /l= o(p) = n( P ), 

e a ,/,=o(p) = e fe)fe= o(p) = e(p). 



(13) 
(14) 

(15) 
(16) 

(17) 

(18a) 
(18b) 
(18c) 
(18d) 
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We calculate in our previous papeJ^S(p), II(p) and the spectrum e(p) up to the third order in J/D. 

T = and h 7^ 0. Taking into account that [H, J2i = an d using Eq. (§ one concludes that h and 
role of chemical potentials for a and b particles, respectively, so that 

G 0i h(w,p) = G(cj + /i,p), 
G 6 , h (w,p) = G(cj - h,p), 
F h {u,p)=F(uj + h,p), 

and, correspondingly, £ a)ft (w,p) = S(w + ft,p), p) = E(w - ft,p), IIh(u>,p) = U(uj + h,p), 

e Q ,fc(p) = e(p) - ft, 
eb,h(p) = e (p) + ft. 



-h play the 

(19a) 
(19b) 
(19c) 



(20a) 
(20b) 



where e(p), £(p), H(p) , G(p ) and F(p) are defined in Eqs. (18 1. 

It is seen from Eq. (20a) that the spectrum of a particles has a gap which vanishes and the spectrum becomes 
unstable at ft < h c \ (T — 0), where 



ft cl (T = 0) = e(p o ) 



(21) 



and po is the momentum at which e(p) has a minimum. This instability signifies a transition to another phase. One 
has po = (71", 7r, 7r) in DTN because exchange couplings arc antiferromagnetic. 

One concludes from Eqs. ( 19 )— ( 20 1 that magnetic field lifts the equivalence between a and b particles. However 
magnetization M(h,T = 0) = (5f)remains zero in the paramagnetic phase as it can be readily seen from Eqs. ([5]), 
(19a I and (19b). Thermal fluctuations make finite the magnetization. 

and h ~ h c i(0). Let us consider the dispersion equation (17) for a particles. We have found its solution in 



our previous paper^l at T = up to the third order in J/D. The aim of the present discussion is to find temperature 
corrections to the spectrum 6T^a,h(p) m the first nonvanishing orders in J/D and T considering the temperature to 
be small enough T -C J. Then, it is convenient to represent the spectrum and self-energy parts using Eqs. ( 19 1 in the 
following form: 



£a,h(p) = £ (P) - ft + 5 T(-a.h{P), 

dZ(u,p) 



£ ,h(ea,fc(p),P) = E ( £ (P)>P) 



*T£ a ,fe(p) + S T ^a,h( e (p) ~ ft' P) 



(22) 
(23) 



e(p) 



where e(p) and p) are the spectrum and the self-energy part at T = 0, 5T^a,^(w, p) is the temperature correction 
to E a ^(w, p) and expressions for e a (p), p) ancl II/i( e a(p),p) can be written similar to Eq. (23). Substituting 

these equations into Eq. (17) and using results of our previous calculation 2 of self-energy parts and e(p) up to the 
third order in J/D we have in the first order in J/D Srfa.p = ^T^a,h{ e (p) — ft, P) an d h c \(T) — e(po) + <5t£ q (0, k ). 
The first order correction in J/D to <5tS q (£(p) — ft, p) is given by the Hartree-Fock diagram shown in Fig. [lja). As 
a result one obtains 



6 T e a<h (p) = 4r a (e(p) - ft, 0, p)M(h, T), 
ft cl (T) = e(p ) +4r Q (0,0,k )M(/i,T), 



where r a (w,p,q) is the vertex and 



ft) 



(24) 
(25) 



(26) 



is equal to the magnetization in the second order in J/D at T <C e(0) + ft. It is explained in our previous papei!^ 
that ladder diagrams give the main contribution to the vertex leading to the Bethe-Salpeter equation for r a (o;,p,q) 
shown in Fig. fl[b) . To calculate the vertex in the leading order in J / D one can use Green's function in the form 
G a (uj, p) — l/(u) — £i(p) + h + iS). When lu ~ J, the solution can be tried in the form r a (cj, p, q) = A(uj) + ( J p — 
J p+q )/4 + B z (uj) -fp + ( q _k )/2 + B xy ^( w )^p+( q _kn-i ii- The solution is quite cumbersome and we do not present it here. 
We point out only that r a ~ J when uj ~ J . 



^k )/2- 

We remind also that the value 



«o = 2r a (0,0,k ) 



(27) 



is an effective two-particle interaction which can be found experimentally at small T as a slope of the plot of ft c i (T) 
vs M C (T) =M(h= h cl {T),T) (see Eqs. pi]) and pel). 



a) 



q + p 



b) r a = 



k -P 






FIG. 1. (a) The Hartree-Fock diagram given the first order correction in J/D to the self-energy part, (b) Diagram equation 
for the vertex F a (uj, p, q) which is involved in the Hartree-Fock diagram. Lines are Green's function ( |lla[ | of a particles. Black 
dots are bare vertexes given by Eq. ( 8b I . 
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FIG. 2. (Color online) Plots of (a) h c \ and (b) M c (magnetization at h — h c i(T)) vs r 3 / 2 in DTN. Circles and diamonds are 
experimental data from RefP and RefP, respectively. Lines are drawn using Eqs. ( |25[ ), ( |26[ ) and parameters Q. 



The specific heat can be obtained from the formula C(h,T) = d(H)/dT using Eq. ([I]) with the following result in 
the first order in J/D: 



C(h,T) 



dT 



(i£(ei(k)-fc)jvMk)-/0 



(28) 



where ei(k) and £3(k) are spectra at T = in the first and in the third orders in J/D, respectively. 

Application to DTN. Equations for M(h,T), h c i(T) and C(h,T) obtained above are applicable at h m /i c i(0) only 
at small enough T so that the thermal corrections to be small. In the case of a quasi-lD system it usually means 
that the temperature cannot exceed the value of the exchange constant between spin chains (~ 0.2 K in DTN). At 
such T the quasi-lD system behaves like a 3D one and we can expect the proportionality to T 3 / 2 of h c i(T) — /i cl (0), 
M C (T) = M(h = h c i(T),T) and C(h = h cl (T),T) expected for QCP of 3D BEC universality class. 
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FIG. 3. (Color online) h c \(T) versus M C (T) in DTN. Circles are exp erim ental data from RefP. Diamonds were putted usin g 
experimental data of RefP for h c i(T) and M C (T) calculated with Eq. (26 1 and parameters Q. Lines are drawn using Eqs. (25 1, 
(26 1 and parameters Q. 
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FIG. 4. (Color online) Specific heat in DTN. Experimental data are taken from RefP^. Lines are drawn using Eq. (28 1 and 
parameters H. 



This proportionality was really observed experimentally in DTN. Fig.[2]shows the experimental data for M C (T) and 
h cl (T) obtained in RefP and RefP for 0.5 K < T < 1 K and 1 mK < T < 300 mK. respectively, together with results 
of our calculations with Eqs. (25), (26) and parameters Q. The agreement between the theory and experiment is 
very good at T < 0.3 K. Large temperature fluctuations come into play at greater T which are not taken into account 
in Eqs. (25) and (26). As a result the deviation from experimental data is noticeable at T > 0.5 K for h c i(T) and at 
T > 0.7 K for MjT). 

Although (and quite expectedly) neither h cl (T) nor M C (T) do not depend linearly on T 3 / 2 at T > 0.7 K in DTN 
(see Fig. 0), it is observed experimentally^ that h c \{T) is a linear function of M C (T) at 0.5 K < T < 1 K as is 
demonstrated in Fig. [3j The effective two-particle interaction is extracted from this plot in RefP as it is explained 
above with the result vq ~ 0.61 meV. Most likely, however, that is renormalized by thermal fluctuations in DTN 
at such large T. We plot in Fig. [3] also h c \(T) vs M C {T) using the low-temperature experimental data of RefP for 
h c i(T) and M C (T) computed from Eq. (26) with parameters Q. The effective interaction obtained in this way is 
equal approximately to 0.44 meV that is 28% smaller than the value experimentally found in RefP at T > 0.5 K and 
that is in excellent agreement with the result of our calculation of vq by Eq. ( [27] ) and parameters Q. 

Experimental data of RefP^for C(h,T) at h « h c i(0) are shown in Fig. [4] together with results of our calculations 
with Eq. (28), parameters Q and expression for es(p) found in RefP. A reasonable agreement between the theory 
and experiment is seen at T < 0.3 K. The specific heat of the model |l]) with parameters Q is proportional to T 3 / 2 
at h = h C 2(T) as well because the spectrum in the fully polarized phase given exactly at T = by e(p) = h — D + J p is 
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also quadratic near its minimum. It is demonstrated in Ref that due to the strong renormalization of the spectrum 
in the paramagnetic phase C(h — h C 2{T),T) is about 6 times larger than C(h = h c i(T),T). Simple calculation of the 
specific heat with the spectrum e(p) shows that in agreement with the experiment C(h — h c2 (T), T) is approximately 
5.7 times larger than Eq. (28) at h = h c i(T) at small T. 
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